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ABSTRACT 



We present a new catalogue of mid-IR sources using the AKARI NEP-Deep survey The InfraRed Camera (IRC) onboard AKARI has 
a comprehensive mid-IR wavelength coverage with 9 photometric bands at 2 - 24 /um. We utilized all of these bands to cover a nearly 
circular area adjacent to the North Ecliptic Pole (NEP). We designed the catalogue to include most of sources detected in 7, 9, 11, 
15 and 18//m bands, and found 7284 sources in a 0.67 deg 2 area. From our simulations, we estimate that the catalogue is ~ 80 per 
cent complete to 200yuJy at 15 - 18//m, and ~10 per cent of sources are missed, owing to source blending. Star-galaxy separation is 
conducted using only AKARI photometry, as a result of which 10 per cent of catalogued sources are found to be stars. The number 
counts at 1 1, 15, 18, and 24 iim are presented for both stars and galaxies. A drastic increase in the source density is found in between 
11 and 15 /um at the flux level of ~ 300yL/Jy. This is likely due to the redshifted PAH emission at 8 fim, given our rough estimate of 
redshifts from an AKARI colour-colour plot. Along with the mid-IR source catalogue, we present optical-NIR photometry for sources 
falling inside a Subaru/Sprime-cam image covering part of the AKARI NEP-Deep field, which is deep enough to detect most of AKARI 
mid-IR sources, and useful to study optical characteristics of a complete mid-IR source sample. 

Key words, infrared: galaxies - surveys - catalogues - methods: data analysis 



1. Introduction 

Over the last two decades, the space infrared missions have 
made deep extragalactic surveys, which produced a growing 
sample of infrared-luminous galaxies. When IRAS discovered 
, a p opulation of ultra-lumin ous infrared galaxies (ULIRGs - 
see ISanders~& Mirabell ll996l for a review), they were regarded 
as very rare objects like QSOs. IRAS also provided good 
evidence for number density evolution of infrared-luminous 
galaxies, although the re dshift range wa s limited to z ^ 0.2 
(Hacking & Houc klll987l: lLonsdale et~ai1ll990t IS aunders et al.1 
1 19901: iGregorich et al.l 119951: iBertin et al.l 119971) . The next in- 
frared mission, ISO, provided further evidence of evolu 
tion towards 



(Elbazetal. 2002; Serjeant et al. 2004; iLe Floc'h et al. 2005; 



Perez- Gonzale z et al. 2005b Babbedge et al. 2006; Caputietal 



1 (iFlores et al. 1999; Oliver et al 



..... . . ^JL 

iRowan-Robinson et al.l 119971: IXul l2000fc ISerieant et al 





1997 




2004 



lOyabu et al.12 005). These measurements are mostly summarized 
in terms of the cosmic star formation rate density, estimated from 
the integrated infrared luminosity function. Probing deeper and 
wider areas, more advanced space telescopes, such as Spitzer, 
AKARI and Herschel, firmly established strong evolution of in- 
frared galaxies, now in the form of evolving infrared luminos- 
ity function. Results from various extragalactic survey fields 
converge specifically in the redshift range z ^ 1. That is, the 
comoving infrared energy density evolves as rapid as (1 + z) 4 



2007; Magnelli etal.1 120 09: Goto et al.l l201d: iRodighiero et al 
120101: iGruppioni et al.l 120101: iMagnelli et al.l 1201 ll) . At z. ~ 1 , 
luminous infrared galaxies (LIRGs) with L IR > 10 n L o are re- 
sponsible for more th an half of this infrared energy density (e.g. 
ILe Floc'h et al.l2005l) . This also indicates that infrared-luminous 
galaxies produced a significant fraction of stellar mass in the 
present universe. Thus, they were indeed major galaxy popula- 
tion in the past, giving important clues as to how galaxies evolve. 
Statistical samples of (U)LIRGs are vital for studies of galaxy 
formation and evolution. 

AKARI is the first Japanese space m ission dedicated to in- 
frared astronomy (Mura kami et al.ll2007l) . AKARI was launched 
in Feb 2006 with the M-V-8 rocket from the Uchinoura Space 
Centre in Japan, as a second generation all-sky surveyor at in- 
frared wavelengths. The first point- source catalogues from the 
all-sky survey were released in March 2010. Along with this 
all-sky survey, AKARI conducted 5088 pointed observations in 
selected areas of sky, during its liquid helium cold phase. Using 
13 percent of the pointed observation opportunities in the cold 
phase, we conducted an extragalactic survey around the North 
Ecliptic Pole (NEP). This NEP survey is two-tiered, consisting 
of the NEP-Deep and the NEP- Wide survey. A salient charac- 
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teristic of this survey is its comprehensive mid-IR wavelength 
coverage - we used 9 photometric bands to span the wavelength 
range from 2 to 24//m. The scientific advantage of obtaining 
such comprehensive wavelength coverage include a capability 
to reliably distinguish starburst-dominated galaxies from those 
with AGN contributions, to obtaining rest- frame mid-IR fluxes, 
free from the uncertainty of complicated ^-corrections. The 
NEP survey provides a valuable input to the study of infrared- 
luminous galaxies. 

This paper describes the mid-IR source catalogue produced 
with the NEP-Deep survey. The description of the data in th e 
NEP- Wide is given elsewhere dLee et al.ll2Q09l:lJeon et al.ll2Q10h . 
We give a brief summary of observations and data reduction 
in section 2, including those for ancillary observations from 
the ground. Source identification and photometry are presented 
in section 3, including performance checks using simulated 
sources. We describe the generation procedures for the cata- 
logue, including optical identification and star-galaxy separation 
in section 4, and give some discussion in section 5. Our summary 
is given in section 6. 

Use of a concordance cosmolog>Q gives a scale of 
1.8 kpc/arcsec at z = 0.1. Therefore, the majority of the galaxies 
appear point-like at z ^ 0. 1 at the resolution of AKARI, that is 
~ 5". 

2. Observations and data reduction 

2.1. AKARI IRC 

The NEP-Deep survey covers a circular field with an area of 
0.67 square degree^ using the IRC onboard AKARI, with a field- 
of-view of 10' x 10'. Si nce details of obse rvations and data re- 
duction are described in Wada et al . (2008), we give only a brief 
summary here. 

The design of AKARFs Sun shield and attitude determina- 
tion system require that the optical axis of the telescope is al- 
ways kept pointing 90° from the direction to the Sun with a 
tolerance of only ±0.6°. Due to this strong visibility constraint 
in its sun- synchronous orbit, deep surveys are possible only 
close to the Ecliptic Poles. Given the presence of the Large 
Magellanic Cloud near the South Ecliptic Pole, we chose the 
North Ecliptic Pole (NE P) for AKARFs unique deep survey field 
dMatsuhara et al.ll2006h . 

The circular field of the NEP-Deep survey can be divided 
into three regions - the central field, and the inner and outer an- 
nuli as shown in Figure [T] The size of the inner annulus is deter- 
mined by the angular separation between the NIR/MIR-S chan- 
nel and the MIR-L channe0. During a pointed observation of 
the inner annulus, the NIR/MIR-S channels always observe the 
opposite side of the MIR-L channel. The angular separation be- 
tween the centers of the IRC channels is 20'. This separation and 
the field-of-view of IRC result in a central hole ~10' in diameter. 
This central hole defines the central region (RA = 17h56m, Dec 
= 66°37 / ) and was covered with one field-of-view of the IRC. 
Because of the difference in the field shape, i.e. a square field- 
of-view for a circular area, some small gaps exist between the 
central region and the inner annulus. Since AKARI has no capa- 
bility to control its roll angle, we utilized the seasonal variation 
of the position angle to cover the survey field. The radius of the 
outer annulus was set to 20' so that there would be no gap be- 
tween the inner and outer annuli. In total we allocated 23, 63, and 

1 a m = 0.3, Q A = 0.7 and H = 70 km sec" 1 MpcT 1 

2 Area covered with L\8W before any masking. 

3 The NIR and MIR-S channels share the same field-of-view. 




Fig. 1. Coverage map of the NEP-Deep survey at 9 jim, showing 
three regions. Square patterns indicate the field-of-view of the 
MIR-S channel of the IRC. See text for details. 



146 pointed observations for the central region, the inner and the 
outer annuli, respectively. We adopted an observing mode with 
neither filter change nor dithering (i.e. so-called AOT05), which 
is optimum for deep surveys. Since each channel of the IRC has 
3 filters, we need at least 3 pointed observations to cover the 
entire wavelength range of the IRC. On average, 4 pointed ob- 
servations per filter were conducted at a given position. 

The observations were executed from 2006 May to 2007 
August. Unfortunately, the quality of the resulting mid-IR im- 
age depends on the season when the observation was conducted. 
From April to August, the observations towards the NEP were 
affected by stray light from the Earth shine. This degradation is 
severer for longer- wavelength images, and affects the north-east 
part of the field in an area of -0.15 deg 2 . 

We used the standard IRC imaging pipeline version 
2007101^3 for each pointed observation, in which dark subtrac- 
tion, subtraction of scattered light inside the camera, correction 
for detector non-linearity, flat fielding, and correction for dis- 
tortion are performed. For the MIR-S and MIR-L images, we 
also removed the diffuse background by subtracting a median- 
filtered image. Astrometry in all of the near-IR and S 1 and S 9 W 
images was determined by the IRC pipeline using 2MASS ob- 
jects as a reference. For images at longer wavelengths, we used 
the coordinates of bright sources in the image one-band shorter 
in wavelength than the image in question. The individual images 
of a given band with calibrated astrometry were combined to 
produce a final mosaicked image by using a publicly available 
software, Swarrfl The images were combined by taking a me- 
dian value of the corresponding individual images. Overviews 
and close-ups of these final images have been presented in Wada 
et al. (2008). 

Wada et al. (2008) have measured the sky noise limit from 
photometry of random positions containing no sources. The 
5<x detection limits of a point source are estimated to be 
9.6/dy in N2 band, 7.5 /Jy (N3), 5.4/Jy (N4), 49 /Jy (57), 
58 jjy (S 9W), 71 /Jy (S 1 1), 1 17 yuJy (L15), 121 jjy (L18 W) and 
276/Jy (L240. Hereafter we use these values as the nominal de- 
tection limits of the NEP-Deep survey. 



4 See AKARI IRC Data User Manual ver 1.3 

5 Available at http : //www . astromatic . net/ 

6 Corrected for the error in the conversion factor used in Wada et al. 
(2008) 
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2.2. Ancillary data 

For the NEP field, we performed various follo w-up obser- 
vation s from X-ray to radio wavelengths (e.g. IWhite et all 
2010a.b). Optical-NIR photometry of mid-IR sources in this pa- 
per utilizes the following data - Subaru/Suprime-cam (BVRi'z'), 
CFHT/Megacam (u*), KPNO 2.1m/FLAMINGOS (JK S ), which 
cover ~ 40 per cent of the NEP-Deep field, except for the u* 
band. An optical source catalogue covering the whole NEP- 
Deep field with C FHT/Megacam (g'r'r'z') is published in 
Hw ang et"aD (120071) . to which readers interested in optical prop- 
erties of more sources are referred. 



2.2.1. Subaru/Suprime-cam 

We used the Subaru Suprime-cam (S-cam) to obtain deep optical 
images for a single field-of-view, i.e. 34' x 27', covering 38 per 
cent of the NEP-Deep field. Observations in the B, R, /', z! bands 
were conducted in June 2003. Additional V-band observations 
were performed in Oct 2003 and June 2004. The central posi- 
tion of the field (RA=17h55m24s, Dec=+66°37'31" [J2000]) 
and the position angle (PA=0°) were selected to avoid the pres- 
ence of bright stars (V < 10 mag) in the field. 

We obtained integration times of 12960 sec for the B band, 
7421 .2 sec for V band, 7200 sec for R band, 6900 sec for V band, 
and 10080 sec for z! band, with a typical seeing of 0.6" - 1". 
Frames which had seeing worse than 1" were not used in pro- 
cessing the final mosaiced images . The data were reduced in a 
stand ard manner using SDFRED (lYagi et all 2002: Ouchi et al. 
2004). Sources are extracted with SExtractor (Bertin & Arnouts 
19961) using the z'-band image as a detection image. The 5-cr 
limiting magnitudes measured with a 2" aperture were 28.4 mag 
for B band, 28.0 mag for V, 27.4 mag for R, 27.0 mag for i' and 
26.2 mag for z! in the AB magnitude system. 



fore we will give only a brief summary here. We conducted 
FLAMINGOS observations in the J and K s bands to cover the 
S-cam field in the NEP-Deep field. We covered the target area 
with observations at 4 different positions. The total effective area 
is 750arcmin 2 , with stellar image sizes ranging from 1.1" to 
2" in FWHM. Data were reduced in the standard manner us- 
ing the IRAF packages. The deepest regions reach / = 22.2 and 
K s = 21.3 AB magnitude in 5 <x. 



3. Source identification and photometry 

3.1. Source detection in MIR-S and MIR-L 

We combined the 57, S9W, and £11 images with the terapix 
software 'SWarp' to produce a detection image for the MIR-S 
channel, median-coadding three separate images. For the MIR- 
L, we used the average of the L15 and L18W images for source 
detection, since the quality of L24 image was systematically 
lower. In the resulting MIR-S and MIR-L detection images, 
the residual from the sky subtraction was reduced significantly, 
which provided more reliable source catalogs. 

For each detectio n image, we ran SExtractor 
(Bertin & Arnoutsl Il996l) to produce initial source lists of 
mid-IR sources. In order to make these initial source lists as 
complete as possible, we adopted a relatively low threshold 
for source detection, i.e. requiring 5 connected pixels with 
> 1.2cr. This resulted in 6746 MIR-S sources and 6719 MIR-L 
sources. Mainly because of the background fluctuation and the 
low threshold, these may include a number of spurious sources. 
Thus, these sources were regarded as 'candidates', and are 
screened with criteria described below. In the final catalogue, 
the MIR-S and MIR-L catalogues were concatenated where 
duplicated entries are removed. 



2.2.2. CFHT/Megacam (u*) 

We obtained an UV image of the whole NEP-Deep field with 
CFHT/Megacam. Observations in the u* band were taken in 
queue mode spread over 12 nights from April 2007 to September 
2007, resulting in a total of 77 frames. The field of view of 
Megacam was large enough to cover the entire NEP-Deep field 
with the central position of RA 17h55h24s and Dec +66°37'32" 
[ J2000] . With the integration time of 600 sec for each frame (2 
frames have 680 sec), we achieved the total integration time of 
773 min. Most of these data-56 frames out of 77-were taken at 
airmasses less than 1.5 and typical seeing of 1.0". 

The data reduction was carried out in a standard manner, 
using software developed for the large-format Megacam data. 
The pre-processing, including bad-pixel masking, overscan and 
bias subtraction, and flat-fielding, were carried out before the 
data delivery using the pipeline system 'Elixir'. By stacking 
the Elixir-processed images, we produced a final mosaiced im- 
age using a software package provided by TERAPIX0 including 
WeightWatcher, SExtractor, SCAMP and SWarp. The depth of 
the final image was estimated to be 24.6 mag [5 <x; AB] from 
photometry of 30000 random positions with 2" aperture. 

2.2.3. KPNO 2.1 m/FLAMINGOS 

llmai et al.l (I2007I) describe observations with KPNO 2.1m 
FLAMINGOS and its data reduction in detail, and there- 



7 http://terapix.iap.fr 



3.2. Photometry and band-merging 

We followed the procedures of Takagi et al. (2007) for photom- 
etry and band-merging for each MIR channel. Starting from the 
initial source position obtained with SExtractor, we searched 
the centroid position for sources in each IRC image. The re- 
evaluated centroid position could shift from the SExtractor posi- 
tion significantly in case of severe source blending or low signal- 
to-noise ratio. Therefore, we performed photometry only if the 
new centroid is less than 3 arcsec away from the initial position, 
which corresponds to ~2 cr of the relative offset between MIR-S 
and MIR-L sources. 

We performed aperture photometry in each IRC image at the 
centroid positions determined above, adopting a 2-pixel (2.9") 
radius for NIR images and 3- pixel (7.0") radius for MIR-S and 
MIR-L images. We then applied the aperture corrections, esti- 
mated from the average digitized point spread function (PSF) 
created from stacking of bright sources around the NEP. In this 
stacking, we did not apply any weight as a function of fluxes, 
since high weight in bright sources also results in high weight 
for associated sources in the outskirts. This reduced the signal- 
to-noise ratio of the outer part of the PSF, and increased the un- 
certainty of the aperture correction. We show the average PSFs 
and the growth curves in Figure[2]and[3l respectively. From these 
growth curves, we derived the aperture correction factors, which 
are tabulated in Table 1 along with other characteristics of the 
PSFs. The aperture-correction factors derived here are defined as 
the ratio of fluxes with two different aperture sizes, i.e. the small 
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Fig. 2. Average PSFs obtained with stacking of bright sources, displayed with a linear grey scale covering 99 % of the dynamic 
range. Solid curves indicate the area with signal-to-noise ratio of 5 or greater. 




Fig. 3. Growth curves of enclosed fluxes as a function of aperture radius. The vertical axis reports relative enclosed fluxes, nor- 
malized to th e value at a 1 5 -pixe l radius. Dotted and dashed lines indicate adopted apertures for photometry and those for flux 
calibration in lTanabe et"aD (120081) . 

Table 1. Summary of point spread functions of IRC 





N2 


N3 


N4 


S7 


S9W 


Sll 


L15 


L18W 


L24 


Reference wavelength \jim] 


2.4 


3.2 


4.1 


7.0 


9.0 


11.0 


15.0 


18.0 


24.0 


# of postage stamps used a 


168 


114 


43 


45 


28 


24 


56 


66 


43 


FWHM (pixel)* 


2.96 


3.06 


2.84 


2.27 


2.19 


2.10 


2.27 


2.48 


2.92 


Aperture radius (pixel) 


2 


2 


2 


3 


3 


3 


3 


3 


3 


Aperture correction factors c 


1.97 


1.96 


1.95 


1.18 


1.21 


1.25 


1.42 


1.53 


1.70 


Pixel scale 


1.46 x 1.46 




2.34 x 2.34 




2.51 x 2.39 





a The number of stacked postage stamp images for PSF creation 
b FWHM computed using the enclosed flux radial profile 
c Correction factors for fluxes 
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apertures adopted here and the IRC calibration aperture^ (10 
pix els for the NIR ch annel and 7.5 pixels for the MIR channels, 
see lTanabe et al.ll2008h . The background level was estimated in 
the sky annulus between a 15- and 20-pixel radius. 

Although this photometric method is simple, it is relatively 
robust against the effects of source blending, compared to pho- 
tometry with larger apertures adopted in the IRC calibration or 
MAG_AUTO photometry in SExtractor. A similar approach was 
successfully appli ed to photome try of blended sources in pre- 
vious works (e.g. iKo et al.l 12009). In heavily confused images, 
photometry with simultaneous PSF-fitting for multiple sources 
may be a better solution. However, there are arguments against 
PSF fitting in our case. It is difficult to accurately determine the 
PSF in each frame, due to uncertainty in the pointing stability 
and the lens aberration whose direction is fixed in the detec- 
tor array. This hampers rigorous PSF-fitting analysis. Moreover, 
faint sources have only a few significant pixels that are useful for 
fitting. Therefore, we adopted aperture photometry, rather than 
PSF-fitting photometry. 

Figure |4] compares aperture- corrected fluxes to those ob- 
tained with the aperture photometry using the aperture size 
adopted for the IRC calibration, i.e. 10 pixels for NIR and 7.5 
pixels for MIR-S and MIR-L. For bright sources, these fluxes 
are consistent with each other, except for the brightest ones that 
have saturated pixels. This confirms that the aperture correction 
obtained from the average PSFs is reasonable. On the other hand, 
fluxes of fainter objects, but detected with ~ 10 <x, show system- 
atic deviation - fluxes with the calibration aperture are system- 
atically larger. For fainter sources, the number density increases, 
and therefore the probability that neighbouring sources fall in the 
large aperture for IRC calibration increases as well. These asso- 
ciated sources would boost the fluxes with large aperture size as 
evidently seen in Figured 

3.3. Astrometry 

The astrometric calibration of IRC images has already been dis- 
cussed by Wada et al. (2008), and will not be repeated here. 
We adopt the coordinate at the shortest wavelength band de- 
tected as the final IRC coordinate of MIR sources. In Figure 
[3 we compare the IRC coordinates with 2MASS coordinates 
for 1178 sources. The mean and cr of the coordinate shift are 
-0.0014" and 0.29" for Right Ascension, and 0.01 1" and 0.32" 
for Declination, respectively. 

3.4. Simulations with artificial sources 

In order to estimate the photometric uncertainties and complete- 
ness of our catalogue, we made a careful simulation to esti- 
mate this. Artificial sources were randomly added to the final 
NEP-Deep images and photometry was done in the same man- 
ner as for real sources. We consider only point sources here, 
since the source detection and photometry are optimized for such 
sources. We used the PSFs shown in Figure [2] for every artificial 
source, adding 200 artificial sources (less than ~3 per cent of 
true sources) at a time, and repeating this for a total of 40,000 
and 160,000 artificial sources, to estimate the photometric errors 
and completeness, respectively. 

8 In calibrating the IRC photometry, Tanabe et al. ( 2008) compared 
the aperture photometry with the expected total fluxes from spectral 
models of calibration stars. Therefore, the IRC calibration factors are 
supposed to include the aperture correction from the calibration aper- 
ture photometry to the total fluxes. 



3.4.1. Photometric errors 

We performed aperture photometry for artificial sources as was 
done for the real sources, and then applied the aperture correc- 
tions. The centroid positions of the artificial sources were de- 
termined through Gaussian-fitting of the simulated point source 
image, using the input position as an initial guess. Thus, the ef- 
fect of noise on the source position was included in the simu- 
lation. In Table 2, we have tabulated the relative flux errors es- 
timated from the input- to-output flux ratio. We note that the ar- 
tificial sources were added to the images with real sources, and 
therefore the photometric errors reported here have some con- 
tribution from the source blending. Here the photometric errors 
were calculated from the scatter of the output fluxes with the re- 
jection of outliers, which are due to the presence of nearby real 
bright sources. Although the catalogue generation method is op- 
timized for point sources, the catalogue includes nearby galax- 
ies, for which IRC fluxes are not as accurate as for point sources. 
We also caution that near-IR fluxes of the several brightest stars 
are saturated and therefore not accurate. 

3.4.2. Completeness and reliability 

The completeness of the resulting catalogue is estimated again 
with the artificial- source method. To reduce the Poisson error 
of the completeness, we produced a large number of artificial 
sources-over 160,000. Since we combine multi-band images 
for source detection, the completeness depends on the colour 
of sources. We consider two cases for the MIR-S; a) flat spec- 
trum in f v and b) red sources, for which we adopt the following 
flux ratios: f^mlfi^m = 6 and f\\ m lfi m = 10 in Jansky units. 
This corresponds to the reddest colours among the sources de- 
tected in all of MIR-S bands0 The average flux ratios of MIR- 
S-detected sources were found to be {f^ml fi^m) = 1-8 and 
(fiifim/fifim) - 2.6. For the MIR-L channel, only the flat spec- 
trum case was considered, because we only combined 15 and 
18yum images where most of galaxies have similar fluxes. 

We ran SExtractor on each simulated image with 200 artifi- 
cial sources. A source was considered to be recovered if its ex- 
tracted position was within a radius of 3". In some other studies, 
the consistency of recovered fluxes is also taken into account 
as a condition of the source recovery, in order to avoid mis- 
identification with real sources. In our case, this flux check is 
not straightforward, since the detection images we made are not 
calibrated. 

We estimated the effect of source confusion in the complete- 
ness analysis as follows. We add artificial sources in the negative 
detection image, created by inverting the background subtracted 
detection image, and compare the results with the case for the 
normal positive image. In Figure we show the completeness 
thus obtained for both positive and negative images. For both 
MIR-S and MIR-L channels, we find that the completeness for 
the negative image is systematically higher than that for the pos- 
itive image. In the completeness analysis, we could wrongly 
identify real sources as artificial sources we generated. If such 
mis-identifications have significant impact on the completeness 
analysis, we would expect that the completeness found for the 
positive image with real sources would be higher than that for 
the negative image. Since this is not the case, the effect of mis- 

9 Our main concern for the source extraction procedure was the ef- 
fects of combining three MIR-S images taken at the different wave- 
lengths. In this procedure, sources with extreme colours may not be 
detected in some of MIR-S bands. In order to evaluate conservatively 
such effects, we conducted the simulation with the most extreme colour. 
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Fig. 4. Fluxes with aperture correction adopted here versus those using the IRC calibration apertures, i.e. 10 pixels for NIR and 
7.5 pixels for MIR-S and MIR-L. Dashed lines indicate 5 cr sensitivity in Wada et al. (2008). Solid lines correspond to the linear 
relation. 




Fig. 5. Comparison of IRC coordinates with the 2MASS coordinates. In the histogram plot, solid and dashed lines indicate ARA 
and ADec, respectively. 



identification is likely to be negligible. The difference between 
the completeness for the positive image and for the negative im- 
age implies another important effect - source blending. When we 
add artificial sources in the positive image, some of them would 
blend with real sources, which results in the positional shift of 
the extracted source and non-recovery. At the 5 cr sensitivity of 
IRC wide band images, i.e. S9W and L18W, the completeness 
for the negative image is ~ 10 % higher than that for the positive 
image in both MIR-S and MIR-L channels. This indicates that 
~10% of mid-IR sources are missed, because of source blend- 
ing. 



In order to evaluate the colour dependence of the complete- 
ness, we calculate the completeness of red sources as a func- 
tion of the average flux of three MIR-S bands, i.e. (/ v )mirs = 
5.67 • fijjm for the assumed colour, which is close to fg m . The re- 
sulting completeness is plotted in Figured and found to be close 
to the case for the flat spectrum. Thus, we conclude that the com- 
pleteness has weak colour dependence, when the average fluxes 
of objects are considered. Since we consider very red sources in 
this simulation, it might be expected that the completeness of the 
longest wavebands, 1 1 /mi, using the multi-band detection im- 
age, could be worse than the case of using the S 1 1 image as the 
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detection image. The completeness of the sin gle-band source ex - 
traction using the same image is reported by lWada et al.l (120081) . 
They obtained 80.7 yuJy for the 50 % completeness at the £11 
band. On the other hand, the 50 % completeness we obtained for 
red sources is 50.6/Jy at the same band, which is comparabl e 
to or even better than the value reported by IWada et al.l (120081) . 
Thus, on average, we see no sign of degraded completeness for 
red sources. 

We estimate the reliability of the catalogue as described 
below. Using the same method as for real sources, we made 
a source catalogue for the negative detection images for both 
MIR-S and MIR-L as described above. In band-merged cata- 
logues, it is expected that sources detected in fewer bands would 
have lower reliability. In Table 3, we summarize the statistics of 
sources with the multiple band detection for 6 bands in MIR-S 
and MIR-L altogether. Although, from the negative image, we 
obtain some spurious sources with multiple-band detection, it 
turns out that almost all of them are rare cases produced by the 
effect of inverted real sources. Therefore, we do not take these 
spurious sources into account, and conclude that sources with 
multiple-band detection are highly reliable. In the negative im- 
ages, there are 113 spurious sources with single band detection. 
These are all faint sources with the fluxes close to the detection 
limits. On the other hand, in the positive images, we find 2138 
sources with single band detection. Given these numbers, we es- 
timate that ~5 % of sources with single-band detection are fake. 
This corresponds to 1.5 % of sources in the final catalogue. 



4. Catalogue 

4.1. MIR-merged catalogue 

At this point, we had two catalogues, based on observations with 
the MIR-S and MIR-L channels, which we then merged together 
by eliminating duplicated entries. In order to choose the best 
photometric results from the duplicated entries, we checked the 
signal-to-noise ratio at each IRC band, and selected the source 
with the maximum number of > 3 cr detections. Finally, we re- 
moved the sources whose fluxes are lower than 5 cr in all MIR 
bands. We adopt the 5 cr sensitivity from Wada et al. (2008). We 
obtained 7284 mid-IR sources in the final catalogue. 

The effects of source blending are not negligible in the pro- 
cess of catalogue generation, since we assume that sources with 
separations of < 3" are the same. There are multiple optical 
counterparts for 29% (8%) of MIR sources within 3" (2"). The 
ratio of source density indicates that only 3% of optical sources 
in the Subaru/S-cam image are detectable in the MIR, and there- 
fore the fraction of sources with serious blending in the IRC im- 
age should be much lower than the source fraction with multiple 
optical counterparts. 

For sources with moderate separation, inspection by eye can 
easily spot the problem of blending which affects both the source 
detection and photometry. In order to call attention to them, we 
set flags (group ID and the number of group member) to identify 
groups of MIR sources with separation of < 5". The sources 
with group ID should be treated with caution. Some of them 
have a common optical counterpart. 

A small portion of this band-merged catalogue is shown in 
Table 4, which is provided in full in the electric version of this 
paper. 
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Fig. 6. Completeness simulation for the MIR-S and MIR-L chan- 
nels. Solid and dashed lines represent the completeness calcu- 
lated with the positive and negative images, respectively. For 
these cases, we assume that artificial sources have a flat spectrum 
in f v . The dot-dashed line in the upper panel indicates the com- 
pleteness for red sources (see text in detail), where we adopted 
the average of 7, 9, and 11 //m flux for red sources. The verti- 
cal dashed lines indicate 5 cr sensitivity in Wada et al. (2008) at 
S9W and L18W bands for MIR-S and MIR-L, respectively. 

4.2. Optical identification 

We made optical identifications within the Subaru/S-cam field. 
Although the S-cam field covers only 38% of the NEP-Deep, it 
is deep enough to detect almost all of the AKARI MIR-detected 
sources. We found multiple optical counterparts for 29 % of MIR 
sources in our 7?-band S-cam image within 3" radius. In order 
to find the best candidate, we adopt the maximum likelihood 
method (Sutherland & Saunders 1992) for the optical identifi- 
cation. The likelihood ratio is defined as 



q(m)f(x,y) 
n(m) 



(1) 



where q(m) is the infinitesimal probability that a MIR source has 
an optical magnitude of m, f(x, y) is the probability distribution 
function of the positional error assumed to be a two-dimensional 
Gaussian, and n(m) is the surface density of background objects 
with magnitude m. To derive the best estimate of q(m), we de- 
fined a sub-sample of 1100 IRC all-band-detected sources, per- 
formed optical identification with a simple nearest neighbour 
method, and visually checked all of the optical identification. 
We derived q(m) using this sub-sample with visually-confirmed 
optical identification. The resulting q(m) is shown in Figure [7] 
along with n(m) from the Subaru 7?-band image. We assume that 
this distribution is not very different from that for general MIR 
sources detected in the NEP-Deep survey. For f(x,y), we adopt 
the astrometric dispersions described above, i.e. 0.29" for RA 
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Table 2. Summary of simulation for photometric errors and completeness 



Flux 








Photometric errors 








Completeness 


(Jy) 


N2 


N3 


N4 


S7 


S9W 


Sll 


L15 


L18W 


L24 


MIR-S 


MIR-L 


3.0e-06 


b 




0.817 














0.03 


0.03 


6.0e-06 




0.678 


0.395 














0.04 


0.03 


1.2e-05 


0.351 


0.326 


0.284 














0.10 


0.04 


2.4e-05 


0.230 


0.152 


0.139 


0.971 


1.056 










0.34 


0.06 


4.7e-05 


0.151 


0.085 


0.088 


0.600 


0.652 


0.834 


1.674 


1.694 




0.70 


0.19 


9.4e-05 


0.089 


0.047 


0.046 


0.326 


0.347 


0.521 


0.781 


0.897 




0.85 


0.54 


1.9e-04 


0.055 


0.031 


0.024 


0.172 


0.183 


0.264 


0.427 


0.456 


1.009 


0.93 


0.79 


3.8e-04 


0.027 


0.015 


0.011 


0.087 


0.094 


0.143 


0.244 


0.261 


0.504 


0.95 


0.89 


7.5e-04 


0.015 


0.008 


0.006 


0.044 


0.047 


0.070 


0.121 


0.137 


0.255 


0.97 


0.94 


1.5e-03 


0.007 


0.004 


0.003 


0.022 


0.024 


0.034 


0.066 


0.069 


0.134 


0.99 


0.95 


3.0e-03 


0.004 


0.002 


0.002 


0.011 


0.012 


0.017 


0.035 


0.036 


0.068 


0.99 


0.97 


6.0e-03 


0.002 


0.001 


0.001 


0.006 


0.007 


0.009 


0.018 


0.017 


0.035 


0.99 


0.98 


1.2e-02 


0.001 


0.001 


0.000 


0.003 


0.003 


0.004 


0.010 


0.010 


0.019 


0.99 


0.99 


2.4e-02 


0.001 


0.000 


0.000 


0.002 


0.002 


0.002 


0.005 


0.005 


0.010 


0.99 


0.99 


4.7e-02 


0.000 


0.000 


0.000 


0.001 


0.001 


0.001 


0.002 


0.002 


0.005 


1.00 


0.99 


9.4e-02 


0.000 


0.000 


0.000 


0.001 


0.000 


0.001 


0.001 


0.002 


0.003 


1.00 


1.00 



a Completeness in MIR-S and MIR-L detection images for the case of constant fluxes. 

b Photometric errors are given if |1 - R ave \ < 0.2, where R ave is the average of the input-to-output flux ratio 



Table 3. The number of sources with multiple MIR band detections 



# of detected bands 


1 


2 


3 


4 


5 


6 


Total 


# of sources 


2138 


1736 


1286 


756 


705 


669 


7284 


Fraction 


0.29 


0.24 


0.17 


0.10 


0.097 


0.092 


1 



and 0.32" for Dec. We calculated L for objects within the search 
radius of 3" and selected the object with the highest L as the 
best candidate for the optical counterpart. A small portion of the 
resulting catalogue of optical counterparts is presented in Table 
5. 

As expected, the resulting optical ID with likelihood ratio is 
sometimes different from the nearest neighbour. For 915 MIR 
sources with multiple optical counterparts in the S-cam field, 
we found that 234 sources have an optical ID which is not the 
nearest neighbour. For these sources, we visually inspected the 
results, but it was not useful to identify correct counterparts. 
Thus, we have serious problems of optical ID for ~ 10% of MIR 
sources. Also, we spotted several pairs of MIR sources which 
happen to have a common optical counterpart. Therefore, optical 
identification of sources with close neighbours should be treated 
with caution. In the catalogue, these sources are identified with 
the group ID flag. 

Out of 7284 mid-IR sources, we found 3162 sources in the 
S-cam field and optical counterparts for all but 79 sources. Most 
of mid-IR sources with no optical ID were detected with only 
a single mid-IR band, which means they are likely to be unre- 
liable sources. However, we find that 9 sources out of 79 were 
detected in more than 3 mid-IR bands. As a result of this visual 
inspection, we found that their optical counterparts are pairs of 
bright sources or lie close to optically bright sources affected by 
spikes, which hamper the correct identification. Interestingly, we 
find a genuine optically-blank mid-IR source even with the deep 
Subaru/S-cam image, which is shown in Figure [8] The R - LI 5 
colour of this source is 9.5 mag (AB), which is about 3 mag red- 
der than the criterion to choose the redde st infrared source s, such 
as faint dust-obscured galaxies (DOGs - iDey et aT] |2008). 



0.30 




R mag 



Fig. 7. Probability distribution function q(m) (shaded histogram) 
and n(m) (empty histogram). Probability densities estimated 
from these histograms are over-plotted with solid curves. 




Fig. 8. Postage stamp of a mid-IR source (MIRS4404) with no 
obvious optical counterparts in the Subaru/S-cam image. Circles 
indicate the IRC/MIR-S position with 3" radius. 
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4.3. Star-galaxy separation 

We mainly use N2 - N3 and N3 - N4 colours for identification 
of stars, since these near-IR bands are the most sensitive IRC 
bands. In the Vega-based magnitude system, the near-IR colours 
of normal stars are close to zero. In this subsection, all magni- 
tudes are given in the Vega-based magnitudes. We firstly made a 
tentative list of stars, based on a ste llarity measured in the CFHT 
r'-band image dHwang etal.i r2007). In Figure we show the 
N3 - N4 versus N2 - N3 colour-colour plot for star-galaxy sep- 
aration. In this plot, objects with a large stellarity, i.e. stellar ob- 
jects, make a clump around zero-colors, as expected. The aver- 
age colours of these stellar objects are (N2-N3) = 0.049+0.057 
and (AG - N4) = 0.006 ± 0.05fJ3. We draw a circular boundary 
in this colour-colour plot to separate stars and galaxies defined 
with the radius AC(= V A (^2 1 N3) 2 + A(N3 - N4) 2 ) of 0.15, 
where we adopt the mean colours of stellar objects as a centre. 
This classification results in 673 stars, out of 7284 sources. 

A large AC for the colour boundary would increase the com- 
pleteness of stars, but could cause mis-classifications of nearby 
galaxies as stars. Therefore, we additionally use following crite- 
ria; -1 <N2- Sll < 1 and N2 < 13 mag for stars. This colour 
is useful to separate stars without mid-IR excess from nearby 
galaxies. 

Furthermore, some bright stars were saturated at the central 
pixels, specifically in the near-IR bands. This causes yet another 
mis-classification if we use near-IR bands for star-galaxy separa- 
tion. To remedy this effect, we adopt additional criteria for stars, 
i.e. SI < 10 mag and -1 < SI - S 11 < 1. With these criteria, 
we finally obtained 720 stars in total. 

There is a prominent galactic source in the NEP-Deep region 
- a planetary nebula NGC6543. The IRC detected filaments of 
this planetary nebula, which should be flagged. Inside this PN 
area, photometry of distant sources is severely affected by these 
filaments. Therefore, we simply flagged out all of sources within 
the PN area. The central position of a circular mask is adjusted 
to encompass the PN area with the radius of 194". 



5. Discussion 

We calculated the number counts of mid-IR sources at 11, 15, 
18 and 24/mi shown in Figure [TO1 In deriving the counts, we 
first defined the area to be used for each image based on the 
weight map. We found most low- weight areas in the north-east 
part of the field, in which AKARI observations suffered from the 
Earth shine. To reduce this effect, we excluded the field with 
the weight of the lowest 20 %. As a result of this cut, we ended 
up with areas of 0.46, 0.53, 0.51, and 0.46deg 2 for 11, 15, 18, 
and 24 /mi, respectively. The excluded area is 0. 167 deg 2 at most. 
The counts are corrected for the complet eness. We compare ou r 
L24 counts with Spitzer 24 /mi counts (iPapovich et al.1 120041) . 
and find a reasonable agreement. However, even with the com- 
pleteness correction, the counts below the 5 cr detection limit ex- 
hibit a significant under-estimation, compared to the Spitzer' s 
counts. Therefore, in Figure [TO) we only show the counts above 
the 5<x detection limits. For other mid-IR number counts i n 
AKARI bands see IWada et al.l (I2007L l2008h : iLee et al.l (120091) : 
iPearson et al . (2010) and Pearson et al. (201 1, in preparation). 

The wavelength-dependence of number counts from 11-to- 
24 /zm could be explained by the SEDs of galaxies, where PAHs 
and hot dust emission dominate. We find the least number of 
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Fig. 9. Near-IR colour-colour plot in Vega magnitudes for star- 
galaxy separation. Objects with large stellarity (> 0.9998) are in- 
dicated with crosses. A large circle indicates the adopted bound- 
ary for star-galaxy sepation (see text in detail). Zero-colours are 
marked with dashed lines. 



sources at 1 1 /mi, where the prominent PAH emission around 
8 jim is redshifted out of the bandpass for z ^ 0.4. 

A typical redshift of catalogued galaxies can be estimated 
from AKARI colours. Figure [TT] shows the colour-colour plot us- 
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Fig. 10. Number counts of galaxies and stars at 11, 15, 18, 
and 24 pm. Large and small symbols with error bars indi- 
cate galaxy and star counts, respectively. See legend for in- 
dividual bands. For the L24 band, the number of stars is not 
enough to derive statistical count s. Dashed lines with symbols 
indicate Spitze r 24 /mi counts (IPapovich et al.1 120041) . Spitzer 
16jimi counts dTeplitz et al.l 120111) . and WISE 12 /mi counts 
(iJarrett et "all 1201 ll) from top to bottom. We show the number 
counts for flux bins containing at least 30 objects and with flux 
greater than the 5 cr detection limit. All counts are corrected for 
completeness. 
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Fig. 11. Colour-colour plot with AKARI bands, using N2, N3,S1 
and S 1 1 . Dots (black) and small crosses (green) indicate galax- 
ies and stars, respectively, detected at more than 3 cr in all of 4 
bands. Solid line (red) indicates the colour of Sc galaxy template 
as a function of redshift. Dot-dashed (black) and dashed (blue) 
lines are for Sb galaxy an d M82 templa t e, re spectively. The SED 
templates are taken from lPolletta et al.l (120071) . 

ing AKARI bands, SI -S 11 versus N2-N3. Both colours have 
a good dynamic range specifically at z < 1, owing to the 1.6/zm 
bump of the stellar emission and the PAH emission at 8 yum. At 
z < 1, N2 - N3 almost continuously increases with increas- 
ing redshift, because of the 1.6/zm bump. On the other hand, 
£7 - S 1 1 has a maximum at z ~ 0.5, since the redshifted PAH 
8 /zm feature is captured by S 1 1 bands at that redshift. These 
trends explain the arch- shaped distribution of galaxies in this 
colour-colour plot. From this colour-colour plot alone, it is safe 
to conclude that most of AKARI mid-IR sources in the NEP- 
Deep field lie at z < 1 . The reddest SI - S 11 galaxies may have 
the stronge st PAH emission f eatures at z ~ 0.5, and are studied 
in detail by Takag i et al.l (1201 Oh as 4 PAH- selected' galaxies. 

6. Summary 

We have generated a mid-IR source catalogue ba sed on im- 
ages from the AKARI NEP-Deep survey presented by lWada et al.l 
(2008). This catalogue was designed to include most of the 
sources detected at 7, 9, 11, 15, or 18 /mi. In the catalogue we 
report 7284 sources in the nearly circular area of the NEP-Deep, 
covering 0.67 deg 2 . From the simulation of artificial sources, we 
have estimated the photometric errors in each band and the com- 
pleteness in both MIR-S and MIR-L channels. The star-galaxy 
separation is based solely on AKARI photometry. For sources in 
the Subaru/S-cam field, where the optical photometry is deep 
enough to detect most of AKARI mid-IR sources, we performed 
optical identification with the likelihood ratio method. The mid- 
IR number counts we derived show a drastic increase of sources 
between 11 /mi and 15 /mi, owing to the effect of redshifted 
PAH emission. At the flux level of 1 mJy, we found the highest 
source density per flux bin in the 24 /mi band. Based on redshift- 
sensitive colours, N2 - N3 and S 1 - S 1 1 , most of AKARI mid-IR 
sources are found to lie at z £ 1 . 
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Table 4. Mid-IR source catalogue of the AKARI NEP-Deep field 



ID 


Name 


RA 


Dec 


N2 


AG 


N4 


SI 


S9W 


sn 


L15 


L18W 


L24 


Star 


Group ID 






[J2000.0] 


fv 


fv 


fv 


fv 


fv 


fv 


fv 


fv 


fv 








AKARI-NEPD 






Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


Jy 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


MIRS464 


J175651.47+661445.1 


269.21447 


66.24587 


1.42e-05 


1.97e-05 


2.16e-05 


4.83e-06 


3.65e-05 


7.28e-05 


O.OOe+00 


9.27e-05 


1.62e-04 








MIRS465 


J175538.70+661455.6 


268.91126 


66.24878 


2.95e-05 


2.62e-05 


2.06e-05 


7.06e-05 


7.86e-05 


1.48e-04 


0.00e+00 


0.00e+00 


0.00e+00 








MIRS466 


J175724.90+661425.9 


269.35379 


66.24055 


1.01e-03 


5.22e-04 


3.35e-04 


1.89e-04 


1.46e-04 


1.87e-04 


3.28e-04 


2.96e-04 


4.04e-04 


1 





MIRS467 


J175712.36+661442.4 


269.30152 


66.24512 


7.83e-05 


6.28e-05 


4.29e-05 


4.51e-05 


1.00e-04 


1.51e-04 


2.74e-04 


3.24e-04 


3.10e-04 








MIRS468 


J175457.95+661440.2 


268.74146 


66.24452 


7.77e-04 


4.55e-04 


2.55e-04 


1.55e-04 


1.04e-04 


7.49e-05 


0.00e+00 


6.96e-05 


1.84e-04 


2 





MIRS469 


J175524.63+661452.7 


268.85266 


66.24798 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


5.16e-05 


7.28e-05 


7.34e-05 


0.00e+00 


0.00e+00 








MIRS471 


J175629.97+661436.9 


269.12488 


66.24360 


6.27e-05 


3.95e-05 


2.93e-05 


1.78e-05 


6.86e-05 


1.52e-05 


3.13e-05 


6.52e-05 


-1.00e+00 








MIRS474 


J175801. 17+661506.9 


269.50489 


66.25192 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


1.24e-04 


9.64e-05 


4.29e-05 








MIRS475 


J175513.85+661417.6 


268.80771 


66.23823 


1.69e-02 


9.98e-03 


5.26e-03 


2.30e-03 


1.74e-03 


l.lle-03 


4.97e-04 


3.28e-04 


8.84e-05 


2 





MIRS477 


J175718.83+661427.3 


269.32849 


66.24094 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


4.65e-05 


0.00e+00 


1.24e-05 


1.47e-04 


2.63e-04 
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MIRS478 


J175721.29+661452.5 


269.33871 


66.24794 


1.07e-05 


1.62e-05 


1.25e-05 


1.32e-05 


0.00e+00 


0.00e+00 


2.08e-04 


2.36e-04 


3.59e-04 
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Notes. — See the electronic version of this paper for the complete catalogue. Flux errors are tabulated only in the electronic version. Column (1): AKARI mid-IR source ID. Column (2): Source Name. 
Columns (3) and (4): AKARI/IRC J2000.0 RA and Dec. Columns (5) through (13): IRC flux densities. Undetected sources have a null value, while -1 means the out-of-field. Column (14): Star/galaxy flag 
[0: galaxies, 1: stars selected with near-IR colours, 2: bright stars identified with mid-IR photometry, 3: sources near the planetary nebula NGC6543]. Column (15): Group ID for close neighbours within 5" 
radius. 



Table 5. Photometry from UV to NIR with CFHT/Megacam, Subaru/S-cam and KPN02.1m/FLAMINGOS 



ID 


RA 


Dec 


AO 


u* 


Au* 


B 


AB 


V 


Ay 


R 


AR 


i' 


Ai f 


z' 


Az' 


J 


AJ 


K s 


AK S 


#ofNeighb. 




Optical / [J2000.0] 


arcsec 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


MIRS1596 


268.44467 


66.39353 


0.89 


24.33 


0.07 


22.80 


0.01 


21.48 


0.01 


20.62 


0.01 


19.79 


0.01 


19.25 


0.01 


99.00 


99.00 


99.00 


99.00 


2 


MIRS1597 


268.35163 


66.39347 


0.27 


23.74 


0.05 


22.84 


0.01 


21.65 


0.01 


21.04 


0.01 


20.60 


0.01 


20.17 


0.01 


99.00 


99.00 


99.00 


99.00 




MIRS1598 


269.13741 


66.39712 


0.10 


22.78 


0.01 


21.97 


0.01 


20.84 


0.01 


20.23 


0.01 


19.87 


0.01 


19.54 


0.01 


18.94 


0.03 


18.46 


0.03 




MIRS1599 


268.94945 


66.39747 


0.06 


99.00 


99.00 


24.23 


0.01 


23.13 


0.01 


22.05 


0.01 


21.39 


0.01 


20.97 


0.01 


20.12 


0.07 


19.27 


0.06 




MIRS1600 


268.73970 


66.39690 


0.26 


23.25 


0.02 


22.75 


0.01 


22.24 


0.01 


21.49 


0.01 


21.23 


0.01 


20.99 


0.01 


20.72 


0.10 


19.94 


0.09 




MIRS1607 


268.80148 


66.39745 


0.18 


23.74 


0.03 


22.87 


0.01 


21.81 


0.01 


21.17 


0.01 


20.86 


0.01 


20.44 


0.01 


20.10 


0.07 


19.22 


0.07 




MIRS1608 


268.99580 


66.39949 


0.12 


22.63 


0.02 


21.62 


0.01 


20.31 


0.01 


19.65 


0.01 


19.31 


0.01 


18.99 


0.01 


18.41 


0.02 


17.99 


0.02 




MIRS1613 


268.68909 


66.40347 


1.95 


22.87 


0.02 


19.07 


0.01 


18.65 


0.01 


18.41 


0.01 


18.83 


0.01 


18.12 


0.01 


18.08 


0.04 


18.56 


0.03 




MIRS1614 


268.97998 


66.39938 


0.13 


22.17 


0.01 


21.70 


0.01 


20.91 


0.01 


20.64 


0.01 


20.52 


0.01 


20.36 


0.01 


20.25 


0.08 


20.18 


0.14 




MIRS1615 


268.66240 


66.38812 


0.44 


19.23 


0.01 


18.52 


0.01 


17.82 


0.01 


17.36 


0.01 


17.17 


0.01 


16.80 


0.01 


16.20 


0.01 


15.90 


0.01 





Notes. — This table includes a subset of MIR sources which are detected with Subaru/S-cam, covering part of the NEP-Deep field. See the electronic version of this paper for the complete catalogue. 
Column (1): AKARI mid-IR source ID. Columns (2) and (3): Subaru/S-cam J2000.0 RA and Dec. Column (4): Angular separation between mid-IR and optical position. Columns (5) through (20): 
UV-Optical-NIR magnitudes and errors in AB. Undetected sources have entries "99". Minimum errors are assumed to be 0.01. Not corrected for the galactic extinction. Column (19): Number of optical 
neighbours within 3" radius in the S-cam catalogue. 



